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We review recent progress in numerical simulations with dually transformed SU{2) LGT, starting 
with a discussion of explicit dual amplitudes and algorithms for SU{2) pure Yang Mills in D = 
3 and D = 4. In the D — 3 case, we discuss results that validate the dual algorithm against 
conventional simulations. We also review how a local, exact dynamical fermion algorithm can 
naturally be incorporated into the dual framework. We conclude with an outlook for this technique 
and a look at some of the current challenges we've encountered with this method, specifically 
critical slowing down and the sign problem. 
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1. Introduction: Duality on the lattice 

The general notion of duality in lattice gauge theory has been present in some form since 
the inception of the field'. Abelian dualities possess several generic features [25|, including the 



interchange of high and low temperature regimes between the original model and its dual, and a 
new gauge symmetry possessed by the dual model closely related to the original gauge symmetry 
(sometimes referred to as "self-duality"). While the high-low temperature duality is also present in 
non-abelian dual models [pO[], a gauge symmetry of the dual is not generally present in the non- 
abelian case. Consequently, the known techniques for simulating a gauge field are not directly 
applicable to the dual theory, and thus novel techniques are required in the non-abelian case. For a 
few recent examples of simulations in the dual U{1) case, the reader is referred to [|l^, 22, 23, 24, 



28]. 



As we shall see below, there is a framework [ |19| , [20| ] within which the formal construction of 
non-abelian duals is straightforward. However, at the level of practical computations one is still 
faced with both non-trivial amplitudes (particularly in the D = 4 case) and constraints on allowed 
configurations that can make it difficult to find ergodic moves. Thus, as an alternative means of 
simulation, progress in the non-abelian case has been much slower and to date the ultimate feasi- 
bility of this approach remains an open question. Notwithstanding, we believe the results reviewed 
here (some of them quite recent) indicate definite progress towards evaluating the potential for dual 
non-abelian simulations. 

In the following we review some of the history of the subject and briefly describe some recent 
advances in computational methods for non-abelian dual simulations. We begin by describing the 
form of the duality in the non-abelian case for groups SU (N) in any dimension, and a systematic 
method for defining the dual amplitude in terms of representation-theoretic quantities. For the 
group SU (2), these methods have recently led to an ergodic algorithm and the derivation of explicit 
amplitudes in both D = 3 and D = 4. Further, the dual degrees of freedom for the pure gauge field 
(closed branched surfaces) couple to the dual polymer picture of fermions in a geometrically natural 
way in the non-abelian case. On a cautionary note, there are features of the resulting dual models 
that present a challenge for practical computations: the dual amplitude is in general not positive 
definite, and extending existing methods to weak coupling is complicated by long auto-correlation 
times. We shall return to these difficulties toward the end of the paper. 



2. Construction of the non-abelian duals 



One of the first explicit constructions of a dual non-abelian model on the lattice was reported 
in the early 1990's by Anishetty et al. for the D = 3, SU (2) case [|l|, ^. While no simulations were 
reported in this work, the dual model was explicitly given and its local, discrete structure clearly 
exhibited; the model was in principle computable provided an ergodic algorithm could be found 
to sample the dual ensemble. A further theoretical development came in 1995, when Halliday and 
Suranyi described the construction of D-dimensional non-abelian dual models on the lattice in 
which the dual degrees of freedom are local. A few years later, an important series of papers by 

' Early forerunners to the dual models discussed below were the strong coupling expansions used for example in | |ll| l 
and discussed further in fnv and references therein. 
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Hari Dass et al. [Q, ^, 10] described the first reported attempts at a full lattice simulation using the 
dual amplitude presented by Anishetty et al. in [|l], ^. 

Further theoretical developments came with the work of Oeckl and Pfeiffer [|^]. Here, a 
general non-abelian gauge theory is considered and the dual theory is constructed in terms of 
representation-theoretic quantities. The degrees of freedom of the resulting dual theory are mani- 
festly local, with a local amplitude defined in terms of contractions of invariant tensors that define 
face, edge, and vertex amplitudes. Borrowing some terminology from loop quantum gravity, this 



structure is sometimes referred to as a lattice spin foam model. The reader is also referred to [19] 
for a clear pedagogical treatment and a generalization of the framework that includes ^-deformed 
and supersymmetric models. 

For concreteness, we shall review here the explicit form of the duality. Recall that given an 
oriented lattice, the action is discretized by forming a sum over plaquettes, S[g\ = T.pepS{gp), 
where the group element gp is the holonomy around an oriented plaquette p, obtained by forming 
a suitable product of the variables ge and their inverses. The lattice partition function is 

3f = f Y[dgee'^>"'''^^^"\ (2.1) 

eeE 

where E is the set of lattice edges. To pass to the dual model, the amplitude at each plaquette is 
expanded into the characters of G: 

e"-'"^ =ZciZi{8p), (2.2) 

i 

where the sum is over all unitary irreducible representations (irreps) of the group G; Xi the 
character function on G associated with the ith irrep. The next step in moving to a dual model is to 
interchange the order of integration and summation: 



UdsAu i^ciXiigp) = £ / ( n n '^■<p)x^{p) tep) ' (2.3) 

K^eeE ) peP i s:P^ / ■' \eeE / peP 

where the maps P — > ^ are labellings of the plaquettes of the lattice by unitary irreducible rep- 
resentations of G. At this stage, the model can already be viewed as dual (the integrals being 
absorbed into the definition of the amplitude). However, this is of limited use computationally as 
the integrals are non-local, high dimensional integrals with no obvious closed form. The key insight 
made explicit by [20] is that the integrals project onto a basis of invariant tensors (also referred to 
as intertwiners). Due to space limitations, we won't derive the rest of the duality transformation. 
The resulting dual model has the following form: 

s:P^ /•' \eeE ) peP fe^peP eeE veV 

where (for D = 3), the amplitudes can be defined as: 




Mf,v)= \J ^ \/f , AE{f,e) = f,4 : ' ^ ^ /,, Ap{f,jp)=cp, (2.5) 
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and ^ is the set of all labellings of the plaquettes by irreps of G and the lattice edges by intertwiners 
(for any plaquette labelling, these have a discrete basis indexed by a finite number of irreps). The 
arguments to the Ay function are as follows: irrep labelling plaquettes incident to v are assigned to 
edges, and intertwiners associated to incident edges are assigned to vertices of the Ay network. In 
the case of G = SU (2), a powerful diagrammatic calculus exists [ |T^ for the explicit evaluation of 
quantities such as Ay (/, v) and Ae (/, ^) ■ It can further be shown that the non-zero contributions to 
the partition function have the geometric interpretation of closed, branched surfaces colored by the 
irreps of G. 

We mention in passing that while an explicit formula for the D = 3 case has been known 
since ^], an explicit form for the D = 4 vertex amplitude practical for simulations was not known 
until very recently with [^. This amplitude (a function of 48 half-integer spin labels) takes order 
7"* time and makes it possible to investigate the possibility of dual simulations the four-dimensional 
case, ultimately the most significant for physical models. 

3. An ergodic non-abelian dual algorithm SU (2) case 

In ||3|] the author and collaborators constructed an ergodic algorithm for sampling the dual 
configurations. In the case G = SU (2) and D = 3, we were able to obtain high quality data up to a 
j3 of 2.85. The results of our runs (using 8^ and 16'' lattices) were in statistical agreement with data 
obtained by conventional means, which we used as a check on our calculations. The expectation 
value of spin (averaged over the lattice) was used as an observable; a plot of our simulations 
including residuals is shown in Figure [l[ While [^] focused on the D = 3, SU (2) case, the algorithm 




Figure 1: Conventional and dual results for the 8 lattice. 

developed has a straightforward generalization with respect to group and dimension. Currently 
simulations in D = 4 are being tested. 

While our first tests used spin as an observable, it is straightforward to compute Wilson loop 



observables in the dual; in fact, a method very much analogous to a dual U{1) used by Panero [23] 
can be applied to build up the expectation values of Wilson loop observables. In the abelian case. 
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this method eliminates exponential decrease in signal to noise, and should behave the same in the 
non-abelian case 



4. Dynamical fermions in the dual 

Recall that when coupling gauge fields to fermions on the lattice, an exact integration of the 
fermionic Grassman variables yields an effective action for the gauge fields given by the determi- 
nant of the quark matrix. This determinant is a non-local function of the gauge degrees of freedom 
and consequently a considerable challenge to compute with [Q]. This determinant can in turn be 
expanded into contributions associated to closed polymers on the lattice. The polymer expansion 
can be viewed as a duality for pure lattice fermions on the lattice, and indeed has been pursued in 



that context as an alternative computational framework [jlj, [ISp. The usefulness of this method is 
best established in two dimensions (see, for example 26|; in three and four dimensions, the 
presence of an alternating sign in the amplitude presents a challenge). 

One can show that applying the spin foam expansion of the gauge degrees of freedom 



( |2.4| ) to each term in the polymer expansion results in a configuration space of closed branched 
surfaces ending on the one-dimensional closed polymers. Remarkably^, the resulting amplitude 
that combines factors from the spin foam and polymer expansion is local in the spin foam and 
fermion degrees of freedom. The factoring of the amplitude into plaquette, edge, and vertex ampli- 
tudes remains unchanged, the main difference being that the vertex amplitudes at vertices along the 
fermionic loop are modified by the presence of charge running through that vertex. As in the pure 
gauge theory, the modified vertex amplitude can be defined in terms of a spin network evaluation; 
examples of modified vertex amplitudes are shown in Figure Q 





Figure 2: Charged vertex amplitudes arising in dual dynamical fermion model 



5. Challenges and Outlook 

The sign problem — If the expectation value of the sign of the amplitude is close to zero then 
methods for extracting expectation values (i.e. "sign trick", see |j8|]) which rely on dividing by the 
sign of expectation value will exhibit large numerical errors. In the simulations performed in |^], 
the presence of a negative amplitude had to be accounted for, however the expectation value of the 
sign remained close to unity, decreasing gradually towards weaker coupling. As weaker coupling 
values were used, long auto-correlation times (discussed next) made it difficult to assess the effect 

^Note however an analogous construction for the f/(l) case had been made previously in [|l2]. 
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of mixed signs in the amplitude. The immediate question under investigation is what the status of 
the sign problem is in four dimensions, as a function of coupling. Depending on whether the naive 
"sign trick" method is successful at couplings used for practical simulations, new methods may 
need to be developed. 

Auto-correlation at weak coupling — As one approaches weaker coupling, the dominant con- 
tributions transition from the isolated bubbles characteristic of strong coupling, to a highly disor- 
dered phase of surfaces extending throughout the lattice. This suggests that a cluster algorithm 
capable of evolving state on larger scales may be a feasible strategy. Of particular interest in this 
connection is a recent proposal by Wolff [27] to apply a worm algorithm in a dual setting, which 
demonstrates very promising critical behavior for a Z2 dual model. 

Outlook — Current computational efforts of the author and his group are focused on opti- 
mizing the D = 4 algorithm, addressing the sign problem in various dual models, and developing 
cluster methods to improve auto-correlation times at weak coupling. On the theoretical front, we 
are also developing methods for evaluating G = SU (3) models, in which finding explicit amplitudes 
remains an open problem. In the near future, we will be investigating incorporating supersymmetry 
using dual methods, as well as dual approaches to chiral fermion simulations. 

While we remain cautious in light of the remaining questions with regard to the sign problem, 
the dual approach continues to provide an intriguing alternative perspective on the physics of lattice 
gauge theory. 
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